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We discnss properties of statistical QCD relevant in Fermi phase space model analysis 
of strange hadron prodnction experimental data. We argne that the analysis resnlts in¬ 
terpreted nsing established statistical QCD properties are demonstrating formation of the 
color deconhned state of matter in relativistic heavy ion collisions at highest CERN-SPS 
energies and at BNL-RHIC, comprising deconhned matter composed of nearly massless 
qnarks and glnons, in statistical eqnilibrinm. 

1. INTRODUCTION 

More than 20 years ago, several theoretical meetings held at the University Bielefeld 
have laid the fonndation for the rapid analytical (pertnrbative QCD) and nnmerical (lat¬ 
tice QCD) development of Statistical Qnantnm Chromodynamics. We will illnstrate in 
this paper, nsing as an example the strangeness signatnre of qnark-glnon (color) plasma 
(QGP), the relevance of this work in the analysis of experimental data, and the resnlting 
identihcation of qnark-glnon plasma formation in relativistic heavy ion collisions. 

The discovery of QGP, nnlike the discovery of a new particle, reqnires a change in the 
nnderstanding of fnndamental hadronic degrees of freedom. The signatnre of QGP is very 
hard to identify, since the hnal state observed in the laboratory always consists of the same 
particles, irrespective of the transitional presence of the deconhned state. What changes 
is the detailed composition of the observed prodnced particle abnndance. We consider 
here the disappearance of strangeness snppression, accompanied by a changed pattern of 
(enhanced) strange antibaryon prodnction, with reestablished baryon-antibaryon symme¬ 
try of thermal spectra. All these phenomena offer compelling evidence for qnark-glnon 
plasma formation 

To demonstrate the formation of QGP, we determine the physical properties of the 
sonrce of hadronic particles we are observing, which we associate with the properties of 
statistical QCD. We show that the anomalons yields of strange hadrons originate from a 
state of matter with physical properties of a qnark-glnon hreball. We begin with a short 
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M [GeV] T/T^ 

Figure 1. Left: ai^\iJ,) as function of energy scale /i for a variety of initial conditions. Solid 
line: as{Mz) = 0.1182 (see the experimental point, which includes the error bar at /r = Mz] 
dotted lines: sensitivity to variation of the initial condition. 

Right: as{27rT) for Tc = 0.160 GeV. Dashed line: as{Mz) = 0.119; solid line = 0.1181; dot- 
dashed line = 0.1156. Dotted line: approximate 2-loop solution, given in Eq.[I|, with choice 
A = 150 MeV. 


survey of statistical QCD results we employ, and follow this with a summary of results 
on hadron production and their interpretation in terms of the QGP state. 

2. STATISTICAL QCD 

2.1. Running coupling constant of QCD 

For the purpose of quark-gluon plasma studies, we require the strength of QCD in¬ 
teraction, and the magnitude of quark mass as function of energy scale. The simplest 
way to obtain these results is to integrate the first order renormalization group equations 
using the perturbative definition of the /?, 7 functions, given an initial value of as{M) and 
mi{M). For the determination of the coupling constant, it has become common to refer to 
the value of as{Mz = 91.19 GeV). We use as the central value as{Mz) = 0.1182 ± 0.002. 

To obtain the solid line in figure Fig.|I| (left) the four-loop /3-function was used, in the 
MS modified minimum subtraction scheme. The dotted lines demonstrate considerable 
sensitivity to the initial value as{Mz)- If as{Mz) were larger, the evaluation of the cou¬ 
pling strength in the ‘low’ energy domain /i < 1 GeV, of interest in QGP studies, would 
be divergent, see dotted lines in Fig.|l| (left) above the solid line. The criticism of the per¬ 
turbative study of strangeness production derives from the belief that the strength of the 
QCD interaction is larger than is now known: to put it differently, an essential prerequire¬ 
ment for the perturbative theory of strangeness production in QGP, is the relatively small 
experimental value as{Mz) — 0.118, which has been experimentally established in recent 
years. It is interesting to note in figure Fig.|l| (left) that a 20% reduction in as{Mz) leads 
to a ‘good’ as(0.1 GeV) < 1. 

When studying thermal processes in QGP at temperature T, the proposed interaction 
scale is /r = 2'kT ~ 1 GeV T/Tc, for Tc ~ 160MeV. In Fig.|I] (right), the solid line 
corresponds to an exactly computed as with physical quark thresholds, evaluated at the 
thermal scale, and expressed in terms of T/Tc. The range of experimental uncertainty 
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in as(T), due to uncertainty in as{Mz), is delimited by dashed and dot-dashed lines 
bordering the solid line. A popular approximation of 




(h) - 


2 


2bi In L 
L _ 


L ^ ln(/xVA^), 


( 1 ) 


agrees, using the standard value A^^^ = 205 ± 25 MeV with exact solution for /i > 2mb. 
Presence of lighter quark loops leads to a major deviation at scales of interest to us. 
Eq. ([^) is represented in Fig. ^ (right) by the dotted line, it misses the exact result by 
factor 2 for < T < 1.75Tc, the effective range of observables emerging in SPS and 
RHIC experiments. The experimental error in determination of is today considerably 
smaller. 

The high sensitivity of physical observables to which depend on a”,n > 2, makes 
it imperative that we do not rely on this approximation (dotted line in Fig. E (right)). 
Yet a hxed value as = 0.25 (instead of as = 0.5) derived from this approximation is still 
often used in a study of the QGP phase properties, energy loss of parton jets, charmed 
quark thermalization, strange quark thermal production, etc. Such treatment of QCD 
interaction underestimates by as much as factor four the interaction with the QGP phase, 
and thus the speed of these processes. In most cases, this mundane factor matters, and 
we see that accurate evaluation of ag at the appropriate physical scale is required. 


2.2. Quark—gluon plasma pressure and phase boundary 

The partition function (i.e., pressure) of the quark-gluon phase is obtained after we 
combine quarks, gluons and vacuum contributions: 


— lnZQGp = PQGP — 


rif 




^C2(7rT)^ + yC3 + 


( 2 ) 


We have inserted here the appropriate quark (uf = 2) and gluon degeneracy. The inter¬ 
actions between quarks and gluons manifest their presence aside of the vacuum structure 
effect in the three coefficients Q 7 ^ 1 , see |T0 
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While higher order terms in the perturbative expansion have been obtained, they sug¬ 
gest lack of convergence of the expansion scheme formulated today. Given that the lowest 
order terms are consistent with thermal lattice results we do not pursue this issue further. 
Rather, we show for which parameter set the lattice and perturbative statistical QGD 
results agree with each other and we use these results to understand the properties of a 
QGP fireball. To reproduce the lattice results, near to T = Tc, the only choice we can 
make is in value of the bag constant, S = 0.19 GeV/fm^ . 

Dotted lines in Fig.|^ show the QGP phase pressure condition Pqgp —0 for different 
velocities of expansion, employing an extrapolation to finite baryochemical potential JH 


The solid line denotes, in the /Xb-T plane the pressure balance between an hadronic point 
particle gas and a stationary quark-gluon phase. Dashed line shows the phase boundary 
with an hadronic gas made of finite size hadrons. However, the fluid flow motion of 
quarks and gluons expands the domain of deconfinement exercising against the vacuum 
force originating in the collective expansion velocity fy 


12 . 
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Figure 2. Pressure balance in //b-Fc-plane. Dotted(from right to left): P = 0 at expansion 
velocity = 0,1/10,1/6,1/5,1/4 and 1/3. Solid line, the phase boundary between hadron gas 
and QGP, with point hadrons, dashed the phase boundary with finite size hadrons. 


The condition P = Pp — B 
B = Pp + {Pp + ep) , 

j- tz „ 


0, including the effect of motion reads 
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K = 



(4) 


where subscript p refers to particle pressure and energy, and where we introduced the 
geometric factor k which characterizes the angular relation between the surface normal 
vector and flow direction. 

Expansion beyond P —> 0 is in general not possible. A hreball surface region which 
reaches condition Eq. and continues to flow outwards must be torn apart. This is 
a collective instability and the ensuing disintegration of the hreball matter should be 
very rapid. Thus we hnd, that a rapidly evolving hreball which supercools into the 
domain of negative pressure is in general highly unstable, and we expect that a sudden 
transformation (hadronization) into conhned matter can ensue in such a condition. It 
is important to note that the situation we described could only arise since the vacuum 
pressure term is not subject to how and always keeps the same value. Looking at the 
high how velocity curves in Fig. ^ we see that an exploding QGP hreball only contained 
by the vacuum can supercool to T ~ 145 MeV. 


2.3. Thermal mass and QGP properties 

It has been shown that it is also possible to reproduce the lattice results using hne 
tuned thermal masses (see table I in [jl^). In Fig.|^, we show the light quark (solid thick 
line) and gluon (dashed thick line) thermal masses which were htted to the lattice QGP 
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Figure 

quarks, 


T/T, 

3. Thermal masses fitted to reproduce Lattice-QCD results |14|, thick solid line for 
and thick dashed line for gluons. Thin lines, perturbative QCD masses for as(/r = 2itT). 


data. The perturbative thermal masses, definitions slightly modihed compared to 
are for quarks and gluons, 




(mj) ^ = 27rasT‘^ ^) ’ 


( 5 ) 


and are also shown in Fig.|^, thin lines (dashed for gluons) obtained using as(/i = 27rT), 
see Fig. H (right). 

We conclude that the thermal masses required to describe the reduction of the number 
of degrees of freedom for T > 2Tc are just the perturbative QCD result. Importantly, this 
means that thermal masses express, in a different way, the effect of perturbative QCD, 
and thus for T > 2Tc, we have the option to use Eq. (^, or the more complex thermal 
mass approach. 

The discussion of the properties of QGP is a wide subject which is beyond the scope of 
this report, and we conclude showing the entropy density as function of temperature in 
Fig.i The solid line is for the case of equilibrated quark-light glue system, in the limit 
of vanishing chemical potential. We note that initially the entropy rises faster than the 
asymptotic behavior, since the QCD interactions weaken, and there is an increase in the 
effective number of acting quark and gluon degrees of freedom. Thus, the drop in entropy 
density is considerable when plasma cools and approaches hadronization condition. In 
order to preserve the entropy content in the hreball when the system expands, from 
T ~ 300 MeV towards 150 MeV, a volume dilution by a factor nine must occur. The 
‘pure glue’ case (short-dashed line) contains as expected about half of the entropy when 
comparing at equal temperature, the addition of strangeness adds at the level of 10% to the 
entropy content, at a given temperature. However, as strangeness is produced temperature 
of the isolated fireball cools, and thus there entropy production is even smaller. 


3. STRANGENESS PRODUCTION IN COLLISION PROCESSES 
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T [MeV] 


Figure 4. Entropy density in chemically equilibrated QGP at Aq = 1 as function of tempera¬ 
ture. Solid line Uf = 2, long dashed line Uf = 2.5, short dashed line ‘pure glue’ Uf = 0. 


3.1. Remarks about thermal strangeness production 

Strangeness, and more generally heavy flavor quarks, can be produced either in the 
hrst interactions of colliding matter, or in the many ensuing less energetic collisions. 
Strange quark mass rUs is comparable in magnitude to the typical temperatures reached 
in heavy ion interactions, and the numerous ‘soft’ secondary parton collisions dominate 
the production of strangeness, and naturally, of the light u, d flavors. 

At the time strange flavor approaches chemical equilibrium in soft collisions, the back 
reaction is also relevant. The quantum mechanical matrix element driving a two-body 
reaction must be, channel by channel, the same for forward and backward going reactions 
due to time reversal symmetry. The actual rate of reaction differs, since there are usually 
considerable differences in statistical and phase space factors. However, the forward an 
backward reactions will balance when equilibrium particle yields are established. This 
principle of detailed balance can sometimes be used to evaluate one of the two reactions. 

This is not the place to repeat the well known details of the kinetic theory of strangeness 
production, [Q based on the perturbative cross sections 0 , evaluated with running QCD 
However, we note that the use of scale dependent QCD parameters, 
s amounts to a re-summation of many QCD diagrams comprising 
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parameters, 
and TTis, with p oc 

vertex, and self energy corrections. A remaining shortcoming of thermal production eval¬ 
uation is that up to day, there has not been a study of the next to leading order hnal state 
accompanying gluon emission in thermal processes, e.g., gg —> ss -|- g. In direct parton 
induced reactions, this next to leading order effect enhances the production rate by a 
factor K = 1.5-3. This may cause a corresponding increase also in the thermal rate of 
production, and thus a reduction in the thermally computed chemical equilibration time 
of strangeness and charm. The magnitudes (up to 0.4 mb) of both types of reactions con¬ 
sidered, quark fusion and gluon fusion, are similar. Gluons dominate flavor production, 
because they are so much more likely to collide in the appropriate production channel. 

Given the cross section and collision frequency, the relaxation constant Tg of chemical 
equilibration can be evaluated and the result is that it takes less than 2 fm/c to equilibrate 
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strangeness if T > 300 MeV, and somewhat less time, in the event that next to leading 
order effects play an important role. Despite this clear indication of near strangeness 
eqnilibration in SPS and RHIC environment, there is considerable difference in the com- 
pntational resnlts [p!^-p2[|, originating in differing assnmptions abont degree of chemical 
eqnilibration of glnons in the early stage of the QGP formation. The relaxation time is 4 
times larger shonld glnons be only 50% eqnilibrated. 


3.2. Strangeness background 

One conld argne that if strangeness is to be nsed as a diagnostic tool of qnark-glnon 
plasma, we need to nnderstand this background production rate of strange hadrons. In 
that context, we are interested to measure how often, compared to light quark pairs, 
strange quarks are made. 

Wroblewski proposed to consider the strangeness suppression factor 
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IP. = 


2(ss) 


(uu) + (dd) 


( 6 ) 


Only newly made ss, uii and dd quark pairs are counted. If strangeness were to be as 
easily produced as light u, d quarks, we would hnd Ws ^ 1. To obtain the experimental 
value for Wg, a careful study of produced hadron yields is required. We cite below results 
obtained using a semi-theoretical method ||2^, in which numerous particle yields are 
described within the framework of a statistical model. In elementary collisions pp, pp, 
e’''e“, a value Wg — 0.22 is obtained, strangeness is thus relatively strongly suppressed. 
On the other hand, in nuclear A-A’ collisions iP more than doubles compared to p-p 
interactions, considered at the same energy. 

To explain the 2-2.5-fold strangeness yield increase in a kinetic model of particle pro¬ 
duction requires a shift towards strangeness production of all particle formation pro¬ 
cesses. In other words, when modeling the enhanced strangeness yields within a variety 
of approaches, in each model a new reaction mechanism must be introduced that favors 
strangeness over none-strange processes. Even at this relatively elementary level of count¬ 
ing hadron abundances, new physics must be introduced. In a model with the deconhned 
phase this new reaction mechanism is due to the presence of mobile gluons, which are 
most effective in making strange quark pairs. Moreover, as the conditions created in 
the QGP become more extreme with collision energy, e.g., initial temperature exceeding 
substantially the strange quark mass, we expect an increase in ip. 

Generally hadron cascade models tuned to produce enough singly strange hadrons, 
predict wrong abundances of the rarely produced particles such as 5, G. We are not aware 
of any kinetic hadron model with or without ‘new physics’ that is capable to reproduce 
the pattern of rare hadron production, along with enhancement of strangeness and hadron 
multiplicity. Moreover, if rapidity spectra are modeled, usually the transverse momentum 
spectra are incorrect, or vice-versa. It thus seems impossible, in a collision model based on 
conhned hadron interactions, to hnd sufficiently many hadron-hadron collisions to occupy 
by hadrons the large phase space (high p±, high y) hlled by products of nuclear collision. If 
indeed a non-QGP reaction picture exists to explain heavy ion collision data, the current 
situation suggests that some essential reaction mechanism has been overlooked for 20 
years. In short, hadron models need a lot more effort to reach satisfactory agreement 











with the experimental results, even regarding rather simple observable such as hadron 
multiplicities, transverse energy production, strangeness yield. Strangeness production 
in thermal processes in hadronic gas were studied carefully Hi, and the time scales 
involved established. It is impossible to explain the observed pattern of strange hadron 
production in kinetic thermal model based on conhned hadron interactions. 

Some researchers abandon the kinetic, i.e., collisional theory of particle production, and 
focus solely on the experimental fact that the observed hadronic multiplicities are result of 
a pre-established statistical distribution near to thermal and chemical equilibrium, which 
works quite well. However, it was already pointed out 15 years ago that such a result 
can be naturally explained in terms of a dynamic theory of a transient deconhned state 
hadronizing in a coalescence model . A detailed study of the subtle deviations in hadron 
yields from precise statistical equilibrium yields allows to understand the hadronization 
mechanism |^, and therefore ultimately also to explore the properties of the hadronizing 
QGP state. 


4. HADRON FREEZE-OUT 

4.1. Chemical nonequilibrium in hadronization 

An extra reaction step ‘hadronization’ is required to connect the properties of the 
deconhned quark-gluon matter hreball, and the experimental apparatus. In this process, 
the quark and gluon content of the hreball is transferred into ultimately free howing 
hadronic particles. In hadronization, gluons fragment into quarks, and quarks coalesce 
into hadrons. Color ‘freezes’, quark-gluon plasma excess entropy has to hnd a way to 
get away, so any additional production is hindered. It is far from obvious that hadron 
phase space (so called ‘hadronic gas’) be used consistently to describe the physics of 
thermal hadronization, and we establish now the consistency criterion assuming that 
entropy production is small, or even null, in hadronization of entropy rich thermal QGP 
into entropy poor hadron phases space. 

We consider the Gibbs-Duham relation for a unit volume, e-|-P = Ta + , and 

combine it with the instability condition of dynamical expansion, Eq. (|^): 


0 — P|h + (-P|h + c|h)' 






■ ("■ * ^1.) (■ * i?t) 


Using global variables we obtain |]T3 

,,2 


E 

~S 


= (T|h + 5T|h) 1 + 




1 — vl 


oT — Pb ~ C/h 
a S/b 


(7) 

( 8 ) 


For RHIC, we have 5T|h < 0.4 MeV, considering that pb < 40 MeV and S/b > 100; at 
top SPS energy, we have pb — 200-250 MeV and S/b 25-45, and thus, 5T|h ^ 5 MeV. 
The particular usefulness of Eq. (H) comes from the observation that it implies: 

I >r|h. (9) 

h 

This is a near equality since the geometric emissivity factor k is positive and small, 
especially so at RHIC, and v/ < 1/3. 
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The Gibbs-Duham relation implies: E/S + PV/S = T + 6T > T. One would think 
that the PV term is small, since the pressure is small as the lattice calculations suggest, 
e/P ^ 7. However, this effect can be compensated by large volume of hadronization. The 
volume is a directly measured quantity. The HBT results place a very severe constraint 
on the emitter size of the pion source. Another way, to constraint the magnitude of V, is 
to relate it to the total yield of pions. With these constraints, the PV term is negligible. 
In Eq. (H), due to super cooling, P was negative, and contributed the n term on the right 
hand side. 

Neglecting the influence of PV, we obtain just as in Eq. (H): E/S ~ T. In Eq. (||), 
we have a universal hadronization constraint, which should be satisfied in models claim¬ 
ing to describe hadron production in relativistic heavy ion collisions. This constraint is, 
in our experience, difficult if not impossible to satisfy in chemical equilibrium statistical 
hadronization models. The reason is that, in such an approach, a rather high temper¬ 
ature T ~ 175 MeV is required to accommodate the high intrinsic entropy content of 
hadron source, but this does not drive up sufficiently the energy content, since for pions 
in equilibrium E/S < T. In the super saturated (yq ~ 1.6) pion gas E/S > T. 

4.2. Phase space and parameters 

Our approach is in its spirit a generalization of Fermi’s statistical model of hadron 
production [p7| , p8| , in that the yield of hadrons is solely dictated by the study of the 
magnitude of the phase space available. 

The relative number of final state hadronic particles freezing out from, e.g., a thermal 
quark-gluon source, is obtained noting that the fugacity fi of the Tth emitted composite 
hadronic particle containing fc-components is derived from fugacities and phase space 
occupancies 7 ^: 

Ni oc /i = JJ yfcAfc. (10) 

k^i 

In most cases, we study chemical properties of light quarks u, d jointly, though on occasion, 
we will introduce the isospin asymmetry. As seen in Eq. (pT]), we study particle production 
in terms of five statistical parameters T, Aq, Ag, 7 q, 7 s. In addition, to describe the shape 
of spectra, one needs matter flow velocity parameters, these become irrelevant when 
only total particle abundances are studied, obtained integrating all of phase space, or 
equivalently in presence of strong longitudinal flow, when we are looking at a yield per 
unit of rapidity. 

Assuming a QGP source, several of the statistical parameters have natural values: 

1. Ag: The value of strange quark fugacity Ag can be obtained from the requirement that 
strangeness balances, (ug — Ug) = 0, which for a source in which all s, s quarks are 
unbound and have symmetric phase space, implies Ag = 1. However, the Coulomb 
distortion of the strange quark phase space plays an important role in the under¬ 
standing of this constraint for Pb-Pb collisions, leading to the Coulomb-deformed 
value Ag ~ 1 . 1 . 

2. 7 g: The strange quark phase space occupancy 7 g can be computed within the frame¬ 
work of kinetic theory and 7 g ~ 1. Recall that the the difference between the two 
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different types of chemical parameters A* and 7 * is that the phase space occupancy 
factor 7 j regulates the number of pairs of flavor ‘i\ and hence applies in the same 
manner to particles and antiparticles, while fugacity Aj applies only to particles, 
while A“^ is the antiparticle fugacity. 

3. Aq! The light quark fugacity Aq, or equivalently, the baryochemical potential /Xb, 
regulate the baryon density of the fireball and hadron freeze out. This density can 
vary dependent on the energy and size of colliding nuclei, and the value of Aq is not 
easily predicted. However, it turns out that this is the most precisely measurable 
parameter, with everybody obtaining the same model independent answer, for it 
directly enters all very abundant hadrons. Since T differs depending on strategy of 
the analysis, the value of pb is not so well determined and we recommend that Aq 
be cited instead of pb = 3T In Aq. 

4. 7 q: The equilibrium phase space occupancy of light quarks 7 q is expected to sig¬ 

nificantly exceed unity to accommodate the excess entropy content in the plasma 
phase, 7 q < 7 q = ~ 1 . 6 . 

5. Tf! The freeze-out temperature Tf is expected to be within 10% of the Hagedorn 
temperature Th — 160 MeV, which characterized particle production in proton- 
proton reactions. 

6 . Turning now to the flow parameters: The collective expansion velocity Vc is expected 
to remain near to the relativistic sound velocity, Vc < l/v^S, the natural flow speed 
of information in the QGP phase. There is a longitudinal velocity which is needed 
to describe rapidity spectra, and there is a hadronization surface motion, aside of 
many further parameters one may wish to use to model profile of velocity of flowing 
matter. 


The resulting yields of final state hadronic particles are most conveniently characterized 
taking the Laplace transform of the accessible phase space. This approach generates a 
function which, in its mathematical properties, is identical to the partition function. For 
example, for the open strangeness sector, we find (with no flow). 


C 


o-Ei/Ti 




k£i 


oclnZ^G 


It is important to keep in mind that: 


( 11 ) 


a. Eq. (P1T| ) does not require formation of a phase comprising a gas of hadrons, but is 
not inconsistent with such a step in evolution of the matter; Eq. (|r^) describes not 
a partition function, but just a look-alike object arising from the Laplace transform 
of the accessible phase space; 

b. the final particle abundances measured in an experiment are obtained after all un¬ 
stable hadronic resonances ‘j’ are allowed to disintegrate, contributing to the yields 
of stable hadrons; 
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c. in some experimental data, it is important to distinguish the two light quark flavors, 
for example experiments are only sensitive to and not and an average over 
isospin does not occur. 


The unnormalized particle multiplicities arising are obtained differentiating Eq. (pT]) 
with respect to particle fugacity. The relative particle yields are simply given by ratios 
of corresponding chemical factors, weighted with the size of the momentum phase space 
accepted by the experiment. The ratios of strange antibaryons to strange baryons of same 
particle type are, in our approach, simple functions of the quark fugacities. 

4.3. Strange hadrons at SPS 

We expect in sudden hadronization chemical non-equilibrium at hadron freeze-out. For 
strangeness, 7 s ^ 1, has been seen early on in experimental data, Full chemical 

non-equilibrium has been first noted in the study of the S-Au/W/Pb collisions at 200 A 
GeV [0. Fitting the yields of hadrons observed, it has been reported that the statistical 
significance increased when chemical non equilibrium was introduced. The statistical 
significance is derived from the statistical error: 




- RLp) 
(arL,)^ 


( 12 ) 


It is common to normalize the total error by the difference between the number of data 
points and parameters used, the so called ‘dof’ (degrees of freedom) quantity. For systems 
we study, with a few degrees of freedom (typically 5-15), a statistically significant fit 
requires that y^/dof < 1. For just a few ‘dof’, the error should be as small as y^/dof< 0.5. 
The usual requirement y^ —> 1 is only applying for infinitely large ‘dof’. 

Turning to the Pb-Pb system at 158y4 GeV collision energy, we consider particle listed in 
table ^ top section from the experiment WA97, for > 0.7 GeV, within a narrow Ay = 
0.5 central rapidity window. Further below are shown results from the large acceptance 
experiment NA49, extrapolated by the collaboration to full Att phase space coverage. The 
total error y^ for the two result columns is shown at the bottom of this table along with the 
number of data points ‘V’, parameters ‘p’ used, and number of (algebraic) redundancies 
‘r’ connecting the experimental results. For r 7 ^ 0, it is more appropriate to quote the 
total y^, since the statistical relevance condition is more difficult to establish given the 
constraints, but since y^/(V — p — r) < 0.5, we are certain to have a valid description 
of hadron multiplicities. We will return to discuss the yields of G, G at the end of this 
subsection. 

In second last column, the superscript ‘s’ means that Ag is fixed by strangeness balance 
and, superscript ‘ 7 q’, in two last columns, means that 7q = 7q = is fixed to 

maximize the entropy content in the hadronic phase space. The fits presented are obtained 
with latest NA49 experimental results, i.e., have updated h“/ 6 , newly published 0 yield 
P 2 | , and we predict the A/p ratio, b is here the number of baryon participants, and 
h“ = 7 r“ -|- K“ -|- p is the yield of stable negative hadrons which includes pions, kaons and 
antiprotons. We see, comparing the two columns, that strangeness conservation (enforced 
in second last column) is consistent with the experimental data shown, enforcing it does 
not change much the results for particle multiplicities. 
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Table 1 

WA97 (top) and NA49 (bottom) Pb-Pb 158A GeV collision hadron ratios compared with phase 
space fits. 


Ratios 

Ref. 

Exp. Data 

Pb|«-7q 

Pb|T''> 

S/A 


33 


0.099 ± 0.008 

0.096 

0.095 

S/A 


33 


0.203 ± 0.024 

0.197 

0.199 

A/a 


33 


0.124 ± 0.013 

0.123 

0.122 





0.255 ± 0.025 

0.251 

0.255 

K+/K- 


34 


1.80± 0.10 

1.746 

1.771 

K-/7r- 


35 


0.082±0.012 

0.082 

0.080 

K»/6 


36 


0.183 ± 0.027 

0.192 

0.195 

h-/b 


32 


1.97 ± 0.1 

1.786 

1.818 

0/K- 


37 


0.145 ± 0.024 

0.164 

0.163 

A/p 

y 

= 

0 


0.565 

0.568 

p/vr- 

all y 


0.017 

0.016 






1.6 

1.15 


N;p 

r 


9;4;1 

9;5;1 


The six parameters (T, r>c, Aq, As, 7q, 7s) describing the particle abundances are shown 
in the top section of table ^ Since the results of the WA97 experiment are not covering 
the full phase space, there is a reasonably precise value found for one velocity parameter, 
taken to be the spherical surface flow velocity Vc of the hreball hadron source. 

As in S-induced reactions where As = 1 0, now in Pb-induced reactions, a value 
A^^ ~ 1.1 characteristic for a source of freely movable strange quarks with balancing 
strangeness in presence of strong Coulomb potential [0, page 3568], i.e., with Ag = 1, is 
obtained. Since all chemical non equilibrium studies of the Pb-Pb system converge to the 
case of maximum entropy, we have presented the results with fixed 7q = 7q = The 

large values of 7 q > 1, seen in table confirm the need to hadronize the excess entropy 
of the QGP possibly formed. This value is derived from both the specihc negative hadron 
h~/h abundance and from the relative strange hadron yields. 

The hts shown in table ^ satisfy comfortably the constraint that E/S > T discussed 
in subsection 01. One of the interesting quantitative results of this analysis is shown in 


the bottom section of table the yield of strangeness per baryon, s/b ~ 0.7. We see, in 
lower portion of table ^ that near strangeness balance is obtained as result of the £t. 

The most rarely produced hadron is the triply strange O(sss) and r2(sss) which are the 
heaviest stable hadrons, Mq = 1672 MeV. The phase space for O is ten times smaller 
than that for S at the conditions of chemical freeze-out we have obtained, and any non- 
statistical hadronization contribution would be first visible in the kl and kl production 
pattern. For the parameters as in table the k2 yields reported by the experiment WA97 
are under predicted by nearly factor two. This yield excess originates at lowest m±. The 
‘failure’ of a statistical hadronization model to describe kl (and by 30% kl) yields has 
several possible explanations. 
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Table 2 

Upper section; statistical model parameters which best describe the experimental results for 
Pb-Pb data seen in Fig.|^. Bottom section: energy per entropy, antistrangeness, net strangeness 
of the full hadron phase space characterized by these statistical parameters. In column two, we 
fix As by requirement of strangeness conservation, and in this and next column we fix 7q = 7q. 
Superscript * indicates values which are result of a constraint. 



Ph\7'^^ 

Pb|/‘> 

T [MeV] 

151 ± 3 

147.7 ± 5.6 

Vc 

0.55 ± 0.05 

0.52 ± 0.29 

K 

1.617 ± 0.028 

1.624 ± 0.029 

As 

1.10* 

1.094 ± 0.02 

7q 

y" = 


7s/7q 

1.00 ± 0.06 

1.00 ± 0.06 

E/h[GeV] 

4.0 

4.1 

s/b 

0.70 ± 0.05 

0.71 ± 0.05 

F;/5[MeV] 

163 ± 1 

160 ± 1 

{s-s)/b 

0* 

0.04 ± 0.05 


12 and 12 enhancement could be caused by strangeness pre-clustering in the deconhned 
phase IP). This would enhance all multistrange hadrons, but most prominently the phase 
space suppressed 12 and 12 yields. This mechanism would work only if pairing of strange 
quarks would be signihcant near to phase transition. There is the possibility that distilla¬ 
tion of strangeness followed by break up of strangelets which process could contribute to 12 
and 12 production. The decay of disoriented chiral condensates has also been considered, 

H- _ 

In view of these pre and post-dictions of the 12 and 12 anomalous yield, one should 
abstain from introducing these particles into statistical hadronization model hts. We 
note that the early statistical descriptions of 12 and 12 yields have not been sensitive to 
the problems we described i^ , |39[| . In fact, as long as the parameter 7 q is not considered, it 
is not possible to describe the experimental data at the level of precision that would allow 
recognition of the 12 and 12 excess yield within statistical hadronization. For example, a 
chemical equilibrium £t, which includes the 12 and 12 yield has for 18 fitted data points 
with two parameters a x^/dof = 37.8/16 Such a £t is quite unlikely to contain all 
the physics even if its appearance to untrained eye is suggesting a very good description 
of experimental data. 

4.4. Strangeness at RHIC 

In the likely event that the QGP formed at RHIC evolves towards strangeness chemical 
equilibrium abundance, or possibly even exceeds it, we should expect a noticeable over 
occupancy of strangeness as measured in terms of chemical equilibrium final state hadron 
abundance. Because much of the strangeness is in the baryonic degrees of freedom, the 
kaon to pion ratio should appear suppressed, compared to SPS results. A more penetrating 
effect of the hadronization of strangeness rich QGP at RHIC is the formation of strange 
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baryons and antibaryons. This high phase space occupancy is one of the requirements for 
the enhancement of multistrange (anti)baryon production, which is an important hadronic 
signal of QGP phenomena In particular, we hope that hadrons produced in phase 

space with a small probability, such as G, G, will be observed with a yield above these 
expectations, continuing the trend seen at SPS. 

Many results from RHIC ^snn = 130 GeV run are still preliminary and the following 
quantitative discussion is probably not the hnal word in this matter. However, the results 
we hnd are very interesting, and in qualitative agreement with the sudden QGP break 
up reaction picture. The data are mainly obtained at the central rapidity region where, 
due to approximate longitudinal scaling, the effects of flow cancel and we can evaluate 
the full phase space yields in order to obtain particle ratios. We do not explore trivial 
results such as = 1, since the large hadron yield combined with the flow of baryon 

isospin asymmetry towards the fragmentation rapidity region assures us that this result 
will occur to a great precision. We also do not use the results for K*,K* since these 
yields depend on the degree of rescattering of resonance decay products Ol^ . The data 


we use has been reported in conference reports of the STAR collaboration of Summer 
2001, which as available are combined with data of PHENIX, BRAHMS, PHOBOS, for 
more discussion of the data origin, see |^. We assume, in our £t in table |^, that the 
multistrange weak interaction cascading H ^ A, in the STAR result we consider, is cut by 
vertex discrimination and thus we use these yields without weak interaction corrections. 

We hrst look at the last column in table |^, the chemical equilibrium £t. Its large 


originates in the inability to account for multistrange 


Similar results are pre¬ 


sented in Ref. , in an equilibrium £t which does not include multistrange hadrons. The 
equilibrium £t yields E/S = 159 MeV < T = 183 MeV contradicting the conditions we 
discussed in depth in subsection O. On the other hand, the chemical nonequilibrium 
hts come out to be in perfect agreement with data, and are consistent with the QGP 
hadronization picture since E/S = 163 > T = 158 MeV and 7 s, 7 q > 1- The value of 
the hadronization temperature T = 158 MeV is below the central expected equilibrium 
phase transition temperature, this hadronization temperatures at RHIG is consistent with 
sudden breakup of a supercooled QGP hreball. The inclusion of the yields of multistrange 
antibaryons in the RHIG data analysis, along with chemical non-equilibrium, allows to 
discriminate the different reaction scenarios. 

We look next at the strangeness content, s/b = 6 , in table the full QGP phase 
space would have yielded 8.6 strange quark pairs per baryon at Aq = 1.085, as we will 
show below in Fig.|^. Thus we conclude that 7 ^*^^ = 6 / 8.6 = 0.7. With this value, 
and using the fitted value 7 ^*^ = 2.1, we compute 7 ^‘^/ 7 ^‘^^ = 2.1/0.7 = 3. The fact 
that the strangeness phase space in QGP is not fully saturated is, on a second careful 
look, in qualitative agreement with kinetic theory predictions, adjusting for the observed 
RHIG-130 run conditions. 

The value of the thermal energy content, E /6 = 25 GeV, seen in table ^ is also in very 
good agreement with expectations once we allow for the kinetic energy content, associated 
with longitudinal and transverse motion. The energy of each particle is ‘boosted’ with 
the factor cosh ?/||. For v± = c/\/3, we have xl = 1-22. The longitudinal flow range 
is about ±2.3 rapidity units, according to PHOBOS results. To obtain the the energy 
increase due to longitudinal flow, we have to multiply by the average, J cosh|/||/|/|| — 
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Table 3 

Fits of central rapidity hadron ratios for RHIC = 130 GeV run. Top section: experimental 

results, followed by chemical parameters, physical property of the phase space, and the fit error. 
Columns: data, full non-equilibrium fit, nonequilibrium fit constrained by strangeness conserva¬ 
tion and supersaturation of pion phase space, and in the last column, equilibrium fit constrained 
by strangeness conservation, upper index * indicates quantities fixed by these considerations. 



Data 

Fit 

Fit 

s-s = 0 

Fit®'! 

s-s = 0 

p/p 

0.64 ±0.07 

0.637 

0.640 

0.587 

p/h- 


0.068 

0.068 

0.075 

A/A 

0.77 ±0.07 

0.719 

0.718 

0.679 

A/h- 

0.059±0.001 

0.059 

0.059 

0.059 

A/h- 

0.042±0.001 

0.042 

0.042 

0.040 

S/S 

0.83 ±0.08 

0.817 

0.813 

0.790 

S-/A 

0.195±0.015 

0.176 

0.176 

0.130 

S-/A 

0.210±0.015 

0.200 

0.200 

0.152 

K-/K+ 

0.88 ±0.05 

0.896 

0.900 

0.891 

K-/7r- 

0.149±0.020 

0.152 

0.152 

0.145 

Ks/h- 

0.130±0.001 

0.130 

0.130 

0.124 

G/S- 


0.222 

0.223 

0.208 



0.257 

0.256 

0.247 

n/VL 


0.943 

0.934 

0.935 

T 


158± 1 

158± 1 

183± 1 

7q 


1.55±0.01 

1.58±0.08 

1* 

Aq 

1.082±0.010 

1.081±0.006 

1.097±0.006 

7s 


2.09±0.03 

2.1±0.1 

1* 

As 

1.0097±0.015 

1.0114* 

1.011* 

E/h[GeV] 


24.6 

24.7 

21 

s/b 


6.1 

6.2 

4.2 

S/b 


151 

152 

131 

E/S[MeY] 


163 

163 

159 

xVdof 

2.95/(10-5) 

2.96/(10-4) 

73/(10-2) 


sinh(2.3)/2.3 = 2.15, for a total average increase in energy by factor 2.62, which takes the 
full energy content to ~ 65 GeV as expected. 

We now consider what experimental hadron yield results imply about total strangeness 
yield in the RHIC hreball. First, we sum up the yield of strange quarks contained in 
hyperons. We have in singly strange hyperons 1.5 times the yield observed in A, since 
remain unobserved. Also, accounting for the doubly strange which are half of the all 
S, and contain two strange quarks, we have: 



1.5 ■ 0.059 + 2 ■ 2 ■ 0.195 ■ 0.059 = 0.133 . 
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Allowing for the unobserved f2 at the theoretical rate, this number increases to (s)Y/h“ = 
0.14. Repeating the same argument for antihyperons the result is 0.10. s and s content 
in kaons is four times that in Ks and thus the total strangeness yield is 


(s + s) 

h- 


0.76, 


with 32% of this yield contained in hyperons and antihyperons. Up to this point, the 
analysis is based on direct measurements and established particle yields. 

We now estimate the increase in the ‘strangeness suppression’ factor Wg, Eq. (^. Cor¬ 
recting for presence of K“ among negatively charged hadrons, and assuming that all three 
pions are equally abundant, we hnd: 


TT -|- 7r+ -|- 

The total number of pions produced comprises pions arising from resonance decays and 
from projectile and target fragments. Thus, as little as half of the pions are originating 
in the newly made qq pairs. In the RHIC ^snn = 130 GeV run, Wg — 0.6. The increase 
compared to SPS is largely due to strangeness content in hyperons. Considering that 
^QGP ^ = 130 GeV, there is still space for a further strangeness yield rise at 

highest RHIC energy, and we hope and expect IT’s —>• 1 when initial temperatures rises to 
well above strange quark mass for sufficient length of time. 


4.5. Strangeness as signatnre of deconfinement 

The rate of strangeness production in QGP is sensitive to the temperature achieved at 
the time gluons reach chemical and thermal equilibrium. There is considerable uncertainty 
how short the time required to relax strangeness flavor in the thermal gluon medium is. 
Gonsideration of a small strangeness mass found in lattice studies of strange hadrons 
implies rapid strangeness chemical equilibration. There is also the probable acceleration 
of equilibration due to the next to leading order effects (iP-factor). In view of this, we 
now discuss a benchmark yield of strangeness, assuming that the equilibration process 
leads to near chemical equilibrium conditions for hadronizing QGP. Specihcally, the light 
quark abundance in the quark-gluon plasma phase is considered at the equilibrium yield, 
while the strange quark yield is characterized by the QGP-phase space occupancy before 
hadronization, 

We consider the ratio of equilibrium strangeness density, arising in the Boltzmann gas 
limit, to the baryon density in a QGP hreball yields: 

Ph l> I (p,T2 + pyir^) 


To hrst approximation, perturbative thermal QGD corrections cancel in the ratio. 
TUg = 200 MeV and T = 150 MeV, we have: 

b ® InAq(In Aq)3/7r2 ■ 


For 


(14) 


The relative yield s/6 is in the approximation considered, nearly temperature indepen¬ 
dent, which allows to gain considerable understanding of strangeness production. This 
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Figure 5. Strangeness yield per baryon as function of Aq in equilibrated quark-gluon plasma. 


ratio is mainly dependent on the value of Aq. This light quark fugacity pertinent to 
the hnal state hadrons is usually quite well determined and does not vary depending on 
the strategy of data analysis. In the sudden hadronization model, the value in QGP is 
imprinted on the yields of hnal state hadrons. 

We show in Fig.|^, as function of Aq — 1 (variable chosen to enlarge the interesting 
region Aq —> 1), the expected relative yield per baryon originating in the QGP, dehned 
in Eq. (p!4D with 7 ^*^^ = 1. At top SPS energy, we see that the equilibrium yield is at 
1.5 strange pairs per participating baryon (for Aq ~ 1.5-1. 6 ). In p-p collisions at the 
corresponding energy, the yield is below 0.3 strange pair per participant |^. We note 
that the p-p reaction yield is suppressed by factor 2 compared to the canonical yield 
resulting when one attempts to use statistical canonical description to size up the phase 
space. The actual experimental yield, is 2.5 times the yield in p-p reactions, see table ||, 
this is however only half as large as in an equilibrated QGP, see Fig.^. There is need to 
have 7 ^*^^ — 0.5-0.7, in both p-p reaction and the SPS top energy Pb-Pb reactions. The 
explanation of this is that the QGP system did not get to be hot enough for long enough 
time to hll the small p-p and very large Pb-Pb phase spaces. 

At the RHIG 130 GeV run, the value Aq = 1.09 allows to understand many particle 
yields at central rapidity. We see, in Fig.|^, that specihc strangeness yield in a QGP hreball 
at equilibrium is expected to be an order of magnitude greater than currently observed at 
SPS top energy. However, comparing to general hadron multiplicity, only a factor 1.5-2 
further strangeness enhancement can at most be expected at RHIG, the remarkable feature 
of the RHIG situation is that this enhancement is found in the (multistrange) baryon 
abundance. Specially, given the large strangeness per baryon ratio, Fig.||, baryons and 
antibaryons produced at RHIG are mostly strange . We are not aware of any reaction 
model other than QGP formation and hadronization which could produce this type of 
anomaly. 

While the specihc strangeness yield s/h is a clear indicator for the extreme conditions 
reached in heavy ion collisions, perhaps an equally interesting observable is the occu- 
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Figure 6. Strangeness occupancy 7s ratio HG/QGP in sudden hadronization as function of Aq. 
Solid lines 7 q = 1, and short dashed 7 q = 1.6. Thin lines for T = 170 and thick lines T = 150 
MeV, common to both phases. 


pancy of the hadron strangeness phase space, 7 ^^*^. The interesting result to expect is 
an enhancement due to the need to hadronize, into a strangeness poor phase, the QGP 
strangeness excess. To see this, we compare the phase space of strangeness in QGP 
with that of resulting HG. The absolute yields must be the same in both phases. This 
hadronization condition allows to relate the two phase space occupancies in HG and QGP, 
by equating the strangeness content in both phases. Ganceling the common normalization 
factor T^/(27r^), we obtain: 


^QGP^QGP^^ 


m. 




J^QGP ) l^J^QGP ) — 


(- 


rus 




„HG^HG 


7qAc 

As 


■Fk^ 


% 


A^As 


(15) 


where Fk and Fy are the phase spaces of kaons and hyperons, respectively. We have, 
without loss of generality, followed the s carrying hadrons in the HG phase space, and 
we have, in hrst approximation, omitted the contribution of multistrange antibaryons. 
We now use the condition that strangeness is conserved to eliminate Ag, in Eq. ([T5|), and 
obtain. 


7°'’ yiTTTlFyFWT+TFFv) 


( 16 ) 


In sudden hadronization, 1 /HG^^qgp ^ growth of volume is negligible, ~ 

the temperature is maintained across the hadronization front, and the chemical 
occupancy factors in both states of matter accommodate the different magnitude of the 
particle phase space. In this case, the QGP strangeness when ‘squeezed’ into the smaller 
HG phase space results in 7 ^^/ 7 ^*^^ > 1 • We show, in Fig.|^, the enhancement of phase 
space occupancy expected in sudden hadronization of the QGP. The temperature range 
T = 150 MeV (thick lines) and T = 170 MeV (thin lines) spans the range considered 
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today at SPS and RHIC. The value of 7 q is in range of the chemical equilibrium in HG, 
7 q = 1 (solid lines), to the expected excess in sudden hadronization, 7 q = 1.6 (short 
dashed lines). 

We note that, for the top SPS energy range where Aq = 1.5-1. 6 , sudden hadroniza¬ 
tion analysis of data implies T ~ 150 MeV, 7 q ~ 1.6, and the value of 7 s increases across 
hadronization by factor 2.7. Since the yield of strangeness seen at SPS implies 7 ^*^^ — 0.6, 
this in turn implies 7 ^*^ ~ 1.6 ~ 7 ^^, as is indeed found in hadron production analysis 
in the sudden hadronization picture, table |^. Because accidentally 7 ^^/ 7 q'^ — 1, one 
can also model the hadronization at SPS energy in terms of an equilibrium hadronization 
model. The pion enhancement associated with the high entropy phase can be accommo¬ 
dated by use of two temperatures, one for the determination of absolute particle yields, 
and another for determination of the spectral shape. Such an approach has similar number 
of parameters, and comparable predictive power. 

However, the SPS condition, ~ 1, is not present at the RHIC energy range, 

where the hadron phase space occupancy for strangeness is signihcantly larger than for 
light quarks, see table QGP hadronization dynamics should emerge clearly from the 
study of (multi)strange hadron yields at RHIC. Using the yields of (multi) strange baryons 
and antibaryons it is possible to discriminate against trivial equilibrium hadronic gas 
models. Moreover, at RHIC we hnd a very strong strange baryon and antibaryon yield 
which is not accessible in kinetic parton models. 

5. FINAL REMARKS 

The deconhned thermal phase manifests itself through its gluon content, which gen¬ 
erates in thermal collision processes a clear strangeness hngerprint of QGP. The SPS 
strangeness results decisively show interesting new physics, with a signihcant excess of 
strangeness and strange antibaryons, and spectral symmetry between baryons and an¬ 
tibaryons. We see, at SPS and at RHIC, considerable convergence of the hadron produc¬ 
tion around properties of suddenly hadronizing entropy and strangeness rich QGP. We see 
hadronization into pious, at 7q —> 7q = ~ 1 . 6 , which is an effective way to convert 

excess of entropy in the plasma into hadrons. Because of large value of 7 q the strange 
particle signature of QGP hadronization become more extreme and clear at RHIC, as the 
strangeness excess is covered by the general hadron yield excess at SPS. The systematic 
behavior as function of collision energy, and other collision system parameters, will in 
future provide a cross check of our reaction picture. 

We have not advocated that the understanding of QGP can or should be based on a 
comparison of the A-A collision system with the elementary systems such as p-p or even 
e'''-e“. We simply do not know at this time if elementary interaction system is in any 
fundamental way different. It could well be that once there is enough energy deposition, 
the relativistic Maxwell demon which leads to the formation of the statistical high entropy 
QCD state is operating in both the small and the big system, and that the difference in the 
observables arises due to a change in the internal excitation (temperature), and size (and 
thus lifespan). For example, there would be no outward (transverse) flow of matter and 
little strangeness thermal production expected in p-p system, compared to A-A collision 
system. The smaller phase space of the p-p system is easier to equilibrate than the A-A 
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collision system, leading to further confusion. 

It is important to remember that, not only at RHIC, and in near future at LHC, 
QGP can be studied. It is very probable that the onset of deconhnement occurs at quite 
modest energies, perhaps in collisions of 20-40H GeV Pb projectiles with a laboratory 
target. An alternate method of QGP exploration is the study of the energy domain 
near the transition. Formation of QGP phase is an endothermic process with onset of 
entropy and strangeness production, and in experiments near to the condition for phase 
transformation, one should be able to recognize these properties quite clearly: for example, 
by the onset as function of energy and of reaction volume of multistrange antibaryon 
production. The measurement of yields as function of energy, and volume (excitation 
functions) should provide information about the nature of the transition to this new 
phase of matter. 

In closing, we have analyzed available RHIC and high energy SPS results and have 
demonstrated that these can be easily understood as result of sudden hadronization of a 
rapidly expanding QGP hreball. We have discussed characteristic nonequilibrium aspects 
of strangeness yields associated with the QGP hreball explosion. Our analysis shows in 
detail how and where the hadronic equilibrium models fail to describe the experimental 
data. 
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